Sign-reversal of the in-plane resistivity anisotropy in hole-doped iron pnictides 
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The concept of an electronically-driven breaking of the rotational symmetry of a 
crystal [Ij, without involving magnetic order, has found experimental support in sev- 
eral systems, from semiconductor heterostructures [2j and ruthenates [3j, to cuprate [4] 
and iron-pnictide [5j superconductors. In the pnictide BaFe2As2, such an "electronic 
nematic state" appears above the magnetic transition dome, over a temperature range 
that can be controlled by external strain [6, 7J. Here, by measuring the in-plane re- 
sistivity anisotropy, we probe the electronic anisotropy of this material over the entire 
nematic/magnetic dome, w^hose end points coincide w^ith the optimal superconducting 
transition temperatures. Counter-intuitively, we find that, unlike other materials, the 
resistivity anisotropy in BaFe2As2 changes sign across the doping phase diagram, even 
though the signs of the magnetic, nematic, and orthorhombic order parameters are 
kept fixed. This behavior is explained by the Fermi surface reconstruction in the mag- 
netic phase and spin-fiuctuation scattering in the nematic phase. The unique behavior 
of the transport anisotropy unveils that the primary role is played by magnetic scat- 
tering in the normal state transport properties of the iron pnictides, suggesting a close 
connection between magnetism, nematicity, and unconventional superconductivity. 
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Superconductivity in all iron-based superconductors is 
found in the proximity of magnetic order, which is pre- 
ceded or accompanied by a structural transition (see [8l|9] 
for review). In the family of 74Fe2As2 (^=Ca, Sr, Ba, Eu) 
compounds, the magnetic order below T/v is a uniaxial 
stripe type [10 , with in-plane ordering vectors (tt, 0) or 
(0,7r). The structural transition at Tg > Tn preempts 
the magnetic stripe ordering and reduces the tetragonal 
symmetry of the lattice to orthorhombic. It has been 
suggested by several experiments jH fTTHl3] and theories 
p!4HT7] that the tetragonal symmetry breaking is caused 
by magnetic, i.e. electronic, rather than elastic degrees 
of freedom, giving rise to an electronic nematic phase, 
see Fig. [l] Doping of the parent BaFe2As2 with elec- 
trons (by Co substitution of Fe, BaCol22) or holes (by 
K substitution of Ba, BaK122) suppresses both magnetic 
and nematic orders, revealing superconductivity with a 
maximum Tc close to the edges of the magnetic/nematic 
dome. Fig. [T] 

Below T5, the sample splits into structural twin do- 
mains with different directions of the orthorhombic Qq 
and bo axes (a^ > bo) [E]. Applying a small uniaxial 
strain not only stabilizes one particular domain below T5, 
but also extends the nematic phase to temperatures well 
above Tg , due to the coupling between the strain field and 
nematic fluctuations (red area in Fig. fl]) [13l [191 1201 • In 



these strained samples, the nematic state anisotropy re- 
veals itself not only in the small orthorhombic distortion 
ao — bo <^ ao-\-bo, but, mostly, in a pronounced anisotropy 
of the in-plane electrical resistivity g = p^/pa — l [6, 7, 21j. 
Regardless of the microscopic nature of the nematic or- 
der parameter, (/p, symmetry imposes that close to the 
nematic transition temperature g = k, cp^ where z^: is a 
coupling parameter. Indeed, this relation was verifled ex- 
perimentally in the parent compound of BaFe2As2 [E]. 
Unsurprisingly, in the only system where the electronic 
anisotropy was systematically studied across the nematic 
dome, the ruthenate Sr3Ru207 [3 , g was found to retain 
the same sign across the phase diagram. 

In the BaFe2As2 family, most of the studies on its ne- 
maticity have focused on the electron-doped compounds, 
flnding ^ > in both magnetic and nematic phases 
[H [71 [22]. More recently, experiments in underdoped 
BaK122 found a rather small, but still positive, value 
of g [21]. Here we report that for larger K-doping lev- 
els, not explored in the previous studies, the anisotropy 
g changes sign both near T^ (nematic phase) and at very 
low temperatures (magnetically ordered phase). Since 
the sign of (f is flxed by the small applied uniaxial strain, 
the coupling parameter k. must change sign across the 
magnetic/nematic dome. We interpret this sign-change 
as a result of asymmetric changes in the magnetic scatter- 
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FIG. 1. (Color Online) Top panel: Phase diagram of the 
electron- and hole-doped Bal22 system focusing on the under 
- doped regime with superconducting region shown in green. 
The domain of the orthorhombic (antiferromagnetic) phase 
below transition from the tetragonal phase, Ts, is shown in 
yellow. Nematic fluctuations (red region) persist far above Ts 
and affect the in-plane resistivity anisotropy in strained sam- 
ples. The lower panel shows the maximum in-plane anisotropy 
^max as a function of doping. Note the significant asymmetry 
of ^max(x) and its sign change in the hole-doped regime. 



ing and in the magnetically reconstructed band structure 
for electron and hole doping. This unique complexity of 
the nematic dome observed in the iron pnictides suggests 
that spin-scattering dominates transport in the normal 
state, establishing a possible connection to a magneti- 
cally mediated mechanism of superconductivity. 

Figure [2] presents our main experimental finding, the 
sign reversal of the in-plane resistivity anisotropy. For 
BaK122 samples with x < 0.202, the resistivity is larger 
along the shorter orthorhombic axis, pa < P6, whereas 
for samples with x > 0.235, the longer axis has higher 
resistivity, pa > pb- This is clearly seen in the inset 
zooming in on the structural transition region. It is also 
apparent that, due to the applied uniaxial strain, the in- 
plane anisotropy starts to appear well above the struc- 
tural transition Ts of the unstrained sample. 

In the bottom panel of Fig.|2]we plot the temperature- 
dependent anisotropy ratio, ^, for several compositions. 
We use two characteristic features of the g{T) curves to 
further analyze the data. First, we plot in the bottom 
panel of Fig. [l] the maximum in-plane anisotropy ^max 
for different x. In the upper panel of Fig. |3J we plot 
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FIG. 2. (Color Online) Upper panel: Normalized 

temperature-dependent resistivity, p(T)/p(300K), for the 
BaK122 samples with x = 0.202 and x = 0.235. The red 
lines show the resistivity along the tto-axis (pa) while the blue 
lines show pb. The inset zooms at the structural transition, 
where a clear reversal of the anisotropy from pb > pa for 
X = 0.202 to pb < pa for x = 0.235 occurs. Lower panel: tem- 
perature dependence of g for different K-doping levels. The 
arrows mark the magnetic transition temperatures. 



g{T « Tat), i.e. at temperatures immediately below the 
magnetic transition, which coincides with the structural 
transition temperature of the unstrainned BaK122 sam- 
ples. 

An interesting general picture emerges from the anal- 
ysis of the experimental observations of this and previ- 
ous studies OH [21]. On the one hand, for all electron- 
doped BaCol22 and parent compositions, the resistiv- 
ity anisotropy is positive, g{x) > 0. The maximum 
anisotropy, gmaxix)^ shown in the bottom panel of Fig.fl] 
is peaked at some intermediate electron doping, vanish- 
ing near the edge of the magnetic/nematic dome, due 
to the vanishing orthorhombic distortion. On the other 
hand, on the hole-doped side, ^max(^) remains positive 
up to moderate hole-doping, but decreases by more than 
one order of magnitude, from ^max(^ = 0) ~ +0.3 to 
^max(^ = 0.202) ^ +0.02. Even more surprisingly. 
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FIG. 3. Top panel: doping dependence of the in-plane re- 
sistivity anisotropy at O.OSTat in both electron- and hole- 
doped BaFe2As2. The red line is a guide to the eye. The 
nematic/magnetic and superconducting domes are shown in 
the background for reference only, without temperature scale. 
Bottom panel: theoretical calculation of the maximum in- 
plane resistivity anisotropy in the paramagnetic phase nor- 
malized by the residual resistivity, po, as a function of the 
chemical potential (in units of the Fermi energy). The insets 
show schematically the differences in the scattering rates cor- 
responding to different Fermi velocities at the hot spots for 
electron- and hole-doped systems. When the hot-spot Fermi 
velocity has a larger component along the b (a) direction, 
the magnetic scattering rate is larger along the b (a) direc- 
tion. The hot spots are obtained by displacing the electron 
pocket (yellow ellipse) by (7r,0), making it overlap with the 
hole pocket (green circle). The shaded red area denotes the 
predicted sign reversal of the anisotropy. 



^max(^) changes sign and exhibits a minimum at a higher 
hole doping concentration, ^max(^ = 0.235) ^ —0.034, 
before it eventually approaches zero close to the edge of 
the magnetic/nematic dome, near x = 0.3. 

This sign-change happens at all temperatures, even 
near T/v, as shown in the bottom panel of Fig. |2] and 
in the upper panel of Fig. [3| Furthermore, as it can 
be seen in the bottom panel of Fig. [2J the magnitude 
of the temperature-dependent anisotropy ratio, ^(T), is 
maximal close to T/v in the parent and slightly hole- 



doped compositions, but increases monotonically on cool- 
ing for higher hole doping levels. This signals that the 
low-temperature anisotropic reconstruction of the band 
structure due to long-range magnetic order plays a pro- 
gressively more important role near the hole-doped edge 
of the nematic/magnetic dome. 

Naively, the fact that p^ > pa across most of the phase 
diagram seems surprising, because the ferromagnetic di- 
rection appears to be less conducting than the antifer- 
romagnetic one. From the orbital ordering point of view 
[23H25] one expects the opposite effect. However, a closer 
look at the Fermi surface reconstruction [71 122] and the 
magnetic scattering mechanisms [26l [27] suggests that 
the anisotropy sign is decided by quantitative factors that 
depend on the electronic structure, going beyond general 
order parameter arguments. Moreover, the different con- 
tributions to the transport anisotropy do not necessarily 
compete with each other. To sort them out, we focus 
first on the nematic paramagnetic phase, in which scat- 
tering by magnetic fluctuations dominates and the band 
structure is not yet reconstructed by the magnetic order. 

The experimental doping evolution of q{T ^ T/v), 
shown in the upper panel of Fig. [Sj displays a rather 
monotonic behavior, changing sign for sufficient hole- 
doping levels and vanishing at the edges of the nematic 
dome. Such a sign-change was previously predicted by 
the theoretical model of Ref. [26 for the anisotropic 
magnetic scattering in the nematic phase. The minimal 
model of Ref. ^26 , shown in the Supplementary Material 
to be consistent with first-principles calculations, consists 
of a circular hole pocket at the center of the square-lattice 
Brillouin zone and two elliptical electron pockets centered 
at momenta (7r,0) and (0,7r), which coincide with the 
magnetic-stripe ordering vectors. In the nematic phase, 
the amplitude of the fluctuations around these two order- 
ing vectors becomes unequal [17 , breaking the tetragonal 
symmetry of the system and inducing anisotropic scat- 
tering. In the experiment, the applied strain selects the 
(tt, 0) fluctuations. 

The work in Ref. [26] showed that, in the presence of 
impurities and in the low- temperature limit - but still 
in the paramagnetic phase - the scattering rates along 
a and 6, and consequently the sign of the resistivity 
anisotropy is controlled by the Fermi velocities of the hot 
spots connecting the hole and electron pockets, see Fig.[3| 
Roughly, in electron-doped compounds, the hot-spots 
Fermi velocity has a larger component along the h direc- 
tion, giving rise to a larger scattering rate for electrons 
travelling along this direction, implying p5 > pa- As 
the system is doped with holes, the hole pocket expands 
while the electron pocket shrinks, changing the position 
and the Fermi velocity of the hot spots, which eventually 
acquires a larger component along the a direction, im- 
plying pa > pb' To go beyond the low-temperature para- 
magnetic limit and make closer connection to the exper- 
iment, we used the same model of Ref. 26 and obtained 




FIG. 4. (Color Online) Polarized light microscopy of BaK122 
samples with doping levels as shown in the panels. Im- 
ages were taken at temperatures just above the tetragonal- 
to-orthorhombic structural transition (T>Ts, right) and 5K 
(T<Ts, left), the latter showing formation of structural do- 
mains due to four different orientations of ao and bo axes. 
The difference in color of the domains is due to spectral de- 
pendence of bireflectance, depending on angle between the ao 
direction and the polarization plane of incident white light. 
Since bireflectance is proportional to the orthorhombic distor- 
tion, the contrast naturally vanishes in the tetragonal phase 
(right panels). 



numerically the maximum resistivity anisotropy in the 
paramagnetic phase as function of doping. The results, 
shown in the bottom panel of Fig. [3J agree qualitatively 
with the experimental data, capturing both the electron- 
hole asymmetry of g{T ^ Tn) and its sign-change in the 
hole-doped regime. This suggests that magnetic scatter- 
ing plays a central role in the paramagnetic phase. 

Our data in Fig|2] also show that the sign-change in 
the resistivity anisotropy persists down to T = 0, deep 



inside the magnetic state, where magnetic fluctuations 
are weaker. In this regime, transport should be governed 
by anisotropic reconstruction of the band structure. By 
employing first-principles calculations, we obtained the 
Fermi surface in the magnetically ordered state and eval- 
uated g in the relaxation-time approximation (see Sup- 
plementary Material). 

Our calculations show that, in the parent compound, 
the reconstructed Fermi surface contains not only quasi- 
isotropic hole and electron pockets, but also Dirac cones 
whose crossing points are located slightly below the Fermi 
level. These features are in good agreement with quan- 
tum oscillation measurements [22 . Due to the large 
value of their Fermi velocities, these Dirac cones dom- 
inate the transport in the magnetic state, giving rise to 
an anisotropy in the resistivity that is in agreement with 
the experimental observations for x = 0. When holes 
are added to the system, the chemical potential shifts 
down, eventually crossing the Dirac points. Using the 
values extracted by ARPES [28 and quantum oscillation 
measurements [22 , we estimate this Dirac-point cross- 
ing to happen between 0.03 and 0.1 holes/Fe. For this 
doping concentration, the contribution to g coming from 
the Dirac pockets is vanishingly small, making the con- 
tribution from the remaining pockets dominant. The lat- 
ter depends on details that are beyond the accuracy of 
first-principles calculations, and can either yield a small 
positive or negative g. Yet, our first-principles approach 
correctly captures the tendency of a vanishing g inside 
the magnetically ordered phase in the hole-doped sam- 
ples. 
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FIG. 5. (Color Online) Temperature-dependent resistivity 
along the orthorhombic ao-axis, pao for samples of BaK122 
with four different potassium doping levels in the strain-free, 
twinned, state. The curves are shown using normalized plots, 
p{T) / p (300 K) . Arrows indicate the transition temperatures 
and were determined by the maximum in the derivative of the 
resistivity and from polarized optical microscopy. 

In conclusion, we observe the reversal of the in-plane 



resistivity anisotropy between electron- and hole-doped 
BaFe2As2 compounds. We attribute such behavior to the 
spin-fluctuation scattering in the nematic phase and to 
the Dirac-cones contribution in the magnetically ordered 
phase. Both mechanisms complement arily predict a re- 
sistivity anisotropy that is larger in the electron-doped 
side and continuously smaller in the hole-doped side. Our 
analysis shows that the change in the sign of the resis- 
tivity anisotropy can be completely understood in terms 
of the band structure of the pnictides, without invoking 
a change from itinerant to local moment physics, as sug- 
gested recently by Ref. [29]. Besides shedding light on 
the singular character of the nematic phase in the iron 
pnictides, our observation serves as a unique flngerprint 
for the predominance of electronic scattering due to spin 
fluctuations, rather than phonons, in the normal state 
of these superconductors. Therefore, our study provides 
strong support to the concept that anomalous magnetic 
scattering in the proximity of a putative quantum critical 
point is intimately connected to magnetically-mediated 
unconventional superconductivity [30l |3T] . 



SUPPLEMENTARY MATERIAL 
Methods 

Single crystals of BaK122 with a size up to 7x7x1 
mm^ were grown from FeAs flux as described in Ref. ^ . 
Potassium content in the samples was determined using 
electron probe microanalysis with wavelength dispersive 
spectroscopy (WDS), see [33 for details. Samples, had 
typical dimensions of 0.5 mm wide, 2 to 3 mm long and 
0.05 mm thick, and were cut from cleaved slabs along 
the tetragonal [110] direction (which becomes either the 
orthorhombic a^ or bo axis below Tsm)- Optical imaging 
was performed on the samples while mounted on a cold 
flnger in a continuous flow liquid helium cryostat (allow- 
ing for precise temperature control in 5K to 300K range) 
using a Leica DMLM microscope. 

Polarized light microscopy was used for visualization 
of the structural domains [18 and the sample selection 
was based on the clarity of domains in the image. In 
Fig. |4] we show polarized light images of the area be- 
tween potential contacts in BaK122 crystals used in this 
study, with x=0.108 (non-superconducting), x=0.202 
(Tc=10K), x=0.235 (Tc=26K) and x=0.259 (Tc=28K). 
Images were taken on cooling at a temperature slightly 
above Tg and at base temperature. The highest contrast 
is observed when the sample is aligned with the [100] 
tetragonal direction at a 45^ to the polarization direc- 
tion of incident light (parallel and perpendicular to the 
orthorhombic ao in different domains.) The contrast of 
domain images depends on the quality of the surface and 
the homogeneity of the samples. Domains are observed 
for samples with x ~ 0.26, and are no longer observed for 



samples with higher doping level, x=0.34, which also do 
not show any features associated with structural transi- 
tion in p{T). 

Selected samples were mounted for four-probe electri- 
cal resistivity measurements, with contacts made by sol- 
dering 50 jam Ag wires using low-resistance Sn soldered 
contacts [34]. Initial resistivity measurements on each 
sample were carried out using a flexible wire arrange- 
ment with no strain applied to the sample (free standing 
state). The results are plotted in Fig. [s] Since this mea- 
surement is performed in the twinned state of the sample 
and contains contributions from both components of the 
in-plane resistivity in the orthorhombic phase, we call 
it pt. Samples were then mounted on a brass horseshoe 
straining device, and strain was applied through the volt- 
age contact wires by deformation of the horseshoe, see 
Ref. [71 f35l for the details of the procedure. Strain in this 
conflguration is applied along the tetragonal [110] axis, 
which selects the orthorhombic ao axis as a preferable 
direction upon cooling below T^. Very soft current leads 
apply no strain. The strain was incrementally increased, 
and for each increment, temperature dependent resistiv- 
ity measurements were made and the domain structure 
was imaged to determine the completeness of detwinning. 
Our previous studies using x-ray scattering have shown 
that when crystals are strained to the point at which no 
domains are visible in polarized microscopy, about 90 % 
of the whole bulk of the sample represents the domain 
whose orthorhombic ao axis is oriented along the direc- 
tion of the strain and therefore parallel to the current. 
Therefore in the detwinned state we are predominantly 
measuring the ao axis resistivity, pao^ while in the un- 
strained state of the sample roughly equal contributions 
of both resistivity components are measured. Therefore, 
we may calculate the bo axis resistivity, p5o(T)=2pt(T) - 

Pao(T). 




FIG. 6. The Brillouin zone of BaFe2As2 corresponding to the 
two Fe unit cell. 



Contribution of the spin-fluctuation scattering 

To justify the minimal model developed in Ref. [26], 
and used here to calculate the resistivity anisotropy 
shown in Fig. 3 of the main text, we use flrst-principles 



band structure calculations. The first step is to un- 
fold the first-principles Fermi surface from the two-Fe 
unit cell to an effective one-Fe unit cell, where the two 
magnetic ordering vectors are given by Qi = (tt, 0) and 
Q2 = (O^tt). In the Bal22 system, due to a sizeable hy- 
bridization with the As orbitals, there are two ways to un- 
fold the overlapping and hybridized electron bands. One 
consists of constructing Wannier functions from a par- 
ticular energy window (usually near the Fermi surface), 
symmetrizing them according to the primitive tetrago- 
nal group, and recalculating the band structure using 
the obtained Wannier functions [36]. This method pro- 
vides a reasonable agreement of the folded bands with full 
DFT calculations, but this agreement may be achieved 
not through an appropriate As-induced hybridization, 
but through unphysically large one-electron hopping. The 
other method makes use of the actual symmetry element 
that reduces the larger unit cell to a smaller one: a glide 
plane. This method has been discussed in details in 
Ref. [33, and more briefly in Ref. |38^. It is this un- 
folding procedure that is appropriate for our purposes. 

In the upper panel of Fig. [7| we show the DFT- 
calculated Fermi surface cross-sections by the plane kz = 
0. Because of the body-centered symmetry the point 
Z= (0, 0, tt) is equivalent to (27r, 0, 0) (see Fig. JGJfor nota- 
tions). Therefore, the ellipticities of the electron pockets 
with respect to the point Z are opposite to those with 
respect to the point F. In the lower panel of Fig. [7| 
we show the Fermi surface cross-sections in the unfolded 
band from which the Fermi surface of the upper panel 
originates. In the unfolded zone, the ellipticity of the 
electron pockets with respect to Z is the same as with 
respect to F, but opposite as compared to the eccentric- 
ity with respect to M and A (see the notation in Fig.p|. 
After folding, the point M folds upon Z, and point A 
upon F, thus creating an impression that the ellipticity 
changes sign with k^. 

We note that there are two hole pockets at FZ but 
only one at MA, indicating that the dominating scatter- 
ing processes are between the hole pockets at FZ and the 
electron pockets. This shows that the model of Ref. [26 , 
considering a two-dimensional three band model with a 
circular hole pocket at the center F of the unfolded Bril- 
louin zone and two elliptical electron pockets located at 
X = (tt, 0) and Y = (0,7r), is an appropriate starting 
point. Within this model, the conductivity tensor is given 
by: 
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FIG. 7. Upper panel: First-principles calculated Fermi sur- 
face cross-sections for Bal22 in the folded Brillouin zone. 
Lower panels: the Fermi surface obtained from the unfold- 
ing procedure discussed in the text across kz = (left panel) 
and kz — T^ (right panel). 



solving the Boltzmann equation in the presence of scat- 
tering by impurities and spin fluctuations, as explained 
in details in Ref. [26 . To obtain the results of Fig. 3 , we 
followed the approach of Ref. [26 , considering a putative 
magnetic quantum critical point and calculating numer- 
ically the temperature-dependent resistivity anisotropy. 
We used the same parameters as those of Ref. [26 , vary- 
ing the chemical potential to mimic the effects of hole 
doping and electron doping. For each value of the chem- 
ical potential, we obtained the complete temperature- 
dependence of the paramagnetic resistivity anisotropy 
and then extracted the maximum anisotropy, obtaining 
the results of Fig. 3 . 
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Contribution of the Dirac cones: first-principles 
calculations 



where '^^ a ~ de\^^x/dh 



is the component i of the Fermi 

velocity of band A and /^ ^ is the non-equilibrium elec- 
tronic distribution function for a (unit) electric field ap- 
plied along the j direction. The latter can be obtained by 



In order to address the kinematic effect of the 
magnetically-reconstructed Fermi surface anisotropy, we 
calculated the transport anisotropy using the actual DFT 
band structure and the constant relaxation time approx- 



imation. It has been established that GGA calculations, 
even though they overstimate the long-range ordered 
magnetic moment by about a factor of two, provide a 
better agreement with the experimental Fermi surface 
than DMFT calculations, which yield the correct mag- 
netic moment [22]. This is, incidentally, another mani- 
festation that not only the local fluctuations, accounted 
for in DMFT, but also nematic fluctuations, consistent 
with a short range magnetic order of several lattice pa- 
rameters, are responsible for the reduction of the ordered 
moment [16 . With this in mind, we used full potential 
LAPW bands, with full GGA magnetization, including 
spin-orbit interaction, to compute the Fermi surfaces and 
the Fermi velocities. 




FIG. 8. The Fermi surface calculated deep inside the magnet- 
ically ordered phase for parent (left panel) electron- (middle 
panel) and hole - doped (right panel) compositions. Green - 
electron pockets, blue - hole pockets. 



The Fermi surface calculated deep inside the magnet- 
ically ordered phase of the parent compound, shown in 
the left panel of Fig.jS) displays two sets of quasi-isotropic 
pockets: one for the holes (blue) and one for the electrons 
(green). It also displays barely visible threads, corre- 
sponding to the Dirac cones. They are green, indicating 
that in the calculations the Dirac points are located very 
slightly below, but practically at the Fermi level. Shift- 
ing the Fermi level up by 30 meV transforms the Dirac 
pockets into clearly visible electron tubes (middle panel 
in Fig. |8|, while moving it down (right panel in Fig. [8| 
creates hole tubes instead. As pointed out in Ref. [39] , 
the Fermi velocity in the Dirac cones is high, therefore 
whenever they are present, they dominate the transport. 
Near the Dirac point, the contribution from a Dirac cone 
to transport is Nvp ex |£^|, where E is the energy distance 
from the Dirac point. 

Using an extrafine mesh of about 60000 inequivalent 
k-points (44x44x40 divisions in the full Brillouin zone) 
we have calculated the transport function Nvp{E) as 
a function of energy, separately for the hole bands and 
electron bands. One can see very clearly two regimes: 
one, where both regular and Dirac electrons contribute, 
and the other, where the latter contribution disappears. 
Except very near the Dirac point itself, Dirac electrons 
dominate transport. Moreover, contrary to the assump- 
tion in Ref. [39 , the anisotropy comes mostly from the 
Dirac cones, and not from the "round" pockets. As Fig. [9] 
shows, the total anisotropy is vanishingly small when the 



Fermi level is crossing the Dirac cones right at the Dirac 
point, which happens for E ^ Ep in the calculation. 
Away from this chemical potential, the anisotropy grows, 
retaining its positive value pb/pa > 1^ foi" either hole or 
electron doping away from this point. 

In the calculation, the Dirac point is very near the 
Fermi level for the undoped compound. To determine 
where it is in the real material, we can use the ARPES 
data of Ref. [28] , or preferably, the quantum oscillation 
measurements of Ref. [22 . In the latter work, the au- 
thors observe that in order to bring the Dirac cone area 
into a complete agreement with the experiment, they 
need to shift the Fermi level in their LDA calculations 
up by 30 meV. In the former work, it was observed that 
the Dirac points are located 30 meV below the actual 
Fermi level. At the same time, they observed a factor 
of about three in the effective mass renormalization, im- 
plying that, in the unrenormalized band structure, the 
corresponding shift of the Fermi level would be 90 meV. 
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FIG. 9. The hole per Fe filling factor (lower curve) and the 
resistivity ratio, pa/pb (upper curve) obtained using band 
structure calculation, as described in the text. Note that 
the "as-calculated" Dirac points lie precisely on the Fermi 
level, and the resistivity anisotropy is essentially zero. In 
order to match the calculated Fermi surfaces with ARPES 
and quantum-oscillation experiments, respectively [22] |28] , 
one needs to shift the calculated Fermi level of BaFe2As2, 
correspondingly, by 30 meV and 90 meV shifts. The result- 
ing positions are indicated by arrows. After these shifts are 
applied, the minimal anisotropy (pa/pb ~ 1) appear at the 
left of the arrows, that is to say, in the hole doping regime. 
Note that no symmetry demands that pa/pb be less than 1, 
small variation in the calculated band structure could render 
it slightly larger than 1, as observed in the experiment. 

One could actually restrict the integration to partic- 
ular parts of the Brillouin zone and separate the Dirac 
contribution from the rest. After that, one could shift 
the Dirac bands by the amount suggested by the experi- 
ment, and the other bands by their measured shift. This 
is however unnecessary, because the anisotropy is domi- 
nated by the Dirac band, and shifts of the other bands 



do not really matter. Therefore, we calculated the total 
anisotropy as a function of doping, using straight DFT 
band structure, and simply marked (Fig. [9| the position 
of the Fermi level of the parent compound corresponding 
to the aforementioned 30 meV and 90 meV shifts. Note 
that the minimal anisotropy, which corresponds to 0.03- 
0.1 h/Fe depending on the shift, is nearly zero, but this is 
accidental; it could be a small positive or a small negative 
number (as in the experiment). DFT certainly does not 
describe these compounds at the level of accuracy suf- 
ficient to distinguish between these possibilities. Yet, it 
consistently describes the vanishing resistivity anisotropy 
in the magnetically ordered phase as holes are doped into 
the system. 
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